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Resume 

We study the spectral properties of the dilute Wigner random real 
symmetric n x n matrices -ff„,p such that the entries H„^p{i,j) take zero 
value with probability l—p/n. We prove that under rather general condi- 
tions on the probability distribution of -ff„,p(j, j) the semicircle law is valid 
for the dilute Wigner ensemble in the limit n,p ^ oo. In the second part 
of the paper we study the leading term of the correlation function of the 
resolvent Gn,p{z) = {Hn,p — zl)~^ with large enough \\raz\ in the limit 
p,n ^ oo, p = 0{n"), 3/5 < a < 1. We show that this leading term, when 
considered in the local spectral scale, converges to the same limit as that 
of the resolvent correlation function of the Wigner ensemble of random 
matrices. This shows that the moderate dilution of the Wigner ensemble 
does not alter its universality class. 

1 Introduction 

The initial interest in the spectral theory of large random matrices has been 
motivated by the stochastic approach to the descriptions of the energy spec- 
trum of heavy nuclei (see e.g. the collection of early papers [22]). Later random 
matrices of infinitely increasing dimensions have seen numerous applications in 
various branches of theoretical and mathematical physics such as statistical me- 
chanics of disordered spin systems, solid state physics, quantum chaos theory, 
two-dimensional gravity (see monographs and reviews [Sj [fil [lOl E]). In mathe- 
matics, the spectral theory of random matrices has revealed deep Hnks with the 
orthogonal polynomials, integrable systems, representation theory, combinato- 
rics, non-commutative probability theory and other theories [19]. 

The first result of the spectral theory of large random matrices was obtained 
by E. Wigner in the middle of 50th [25] on the eigenvalue distribution of the 
ensemble An of n x n real symmetric matrices of the form 

Anii,j) = -^a{i,j), i,j = l,...,n, (1.1) 
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where {a{i,j), 1 < i < j < n} are independent random variables. E. Wigner [25] 
proved that in the case when random variables a{i,j) have symmetric probability 
distribution with the second moment v'^ and such that all random variables 
a{i,j) have all moments finite, then the eigenvalue counting function 



a„(A,A„) = ltt{Af^ <A}, (1.2) 



n 



where X\ < ... < A„ denote the eigenvalues of A„, weakly converges in 
average as n ^ oo to the limiting function (Tsc(A), with the derivative a-'scW = 
Psc of the semicircle form 

p,,(A) = a' (A) = /^4z;2-A2, if |A| < 2v; 

27rw2 |o, otherwise. 

This limiting distribution (1.4) is referred to as the Wigner distribution and the 
convergence 

fT„(A,A„) cr,c(A) (1.4) 

is known as the semicircle (or Wigner) law. Also the ensemble of random ma- 
trices {An} (1.1) with jointly independent centered random variables having 
the variance is called the Wigner ensemble. Wigner has proved convergence 
(1.4) with the help of the averaged moments of (7„(A, An) determined in natural 
way by the traces of powers of An . 

Another proof of the semicircle law (1.4) can by obtained in frameworks 
of the resolvent approach introduced first in the random matrix theory by 
V. Marchenko and L. Pastur [18]. Moreover, it can be shown that the norma- 
lized trace of the resolvent gn{z) = ^ Tr {A — z)~^ converges to the Stieltjes 
transform 'w(z) of (Tsc(A) under much more relaxed conditions than those of the 
Wigner's original proof [18t I21j. 

The further progress in the studies of the resolvent of random matrices of the 
Wigner ensemble is related with the asymptotic expansions of the covariance 
function 

Cn{zi,Z2) = E{g„(2;i)5„(z2)} - E{.g„(zi)}E{5„(z2)} (1.5) 

that is sometimes referred to as the correlation function of the resolvent. In 
paper [13], it is proved that if arbitrary distributed random variables have 
the fifth moment finite, then the asymptotic expansion of C„(zi, Z2) is given by 

Cn{zi,Z2) = \f{zi,Z2) + o{\), \\mZj\>2v, (1.6) 

where the leading term f{zi,Z2) depends on the limiting Stieltjes transform 
■w{z) and on the moments Ea(i, j)^ = (1 -f 5ij)v'^ and V4 = Ea(z, j)^ ; the form 
of this term is such that in the local scaling limit the following convergence holds 

1 . 7" , f s 1 

This expression coincides with the averaged version of the density-density cova- 
riance function obtained by F. Dyson for the Gaussian Orthogonal Ensemble of 
random matrices [7]. The right-hand side of (1.7) does not depend on the mo- 
ments V21 of random variables and this result supports the universality conjec- 
ture for the local spectral properties of random matrices in the bulk of the 
spectrum. 



lim — /(A- — +iO,A+ — -iO) = — \\\<2v. (1.7) 
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During two last decades, there is a growing interest to certain versions of the 
Wigner ensemble of random matrices named by the dilute random matrices (for 
example, see the review [ij)- The spectral properties of this kind of ensembles 
have been intensively studied numerically and analytically in theoretical physics 
literature (see [8l[9l[23] for the earlier results and [HI [25 for the recent advances 
and references). In particular, it is shown that in the limit of large n and not 
too strong dilution the semicircle law is valid for the dilute Wigner ensembles 
|23| ; also the universal behavior of the density-density correlation function is 
detected on the theoretical physics level of rigour [9]. 

The aim of the present paper is two-fold. First, we prove the analog of the 
statement (1.5) for the dilute Wigner ensembles of random matrices such that 
a.ij belong to fairly wide classes of random variables. To do this, we develop the 
cumulant expansions approach proposed in [13] to study the resolvent of the 
Wigner random matrices. In the second part of the present paper, we use the 
technique developed and prove analogs of relation (1.6) for the dilute Wigner 
random matrices. We show that in certain asymptotic regimes the analogs of 
the universality relation (1.7) are true. This allows one to conclude about the 
universality of the local spectral statistics of dilute Wigner random matrices. 

The outline of this paper is as follows. In Section 2, we define the dilute 
Wigner random matrix ensemble i/„,p and formulate our main results. In Sec- 
tion 3 we prove the semicircle law. In Section 4 we study the corresponding 
correlation function Cn,p{zi, Z2) of the resolvent Gn,p(z) — {Hn.p — zl)^^ we 
show that the variance of the normalized trace of the resolvent Varg„.p(z) is 
bounded by (np)"^ ; also we find the leading terms of Cn,p(zi, Z2). In Section 5 
we prove the auxihary statements used in Section 4. In Section 6 we study the 
asymptotic properties of the leading terms of Cn,p{zi, Z2) and prove analogs of 
relation (1.7). 

2 Main results and the scheme of the proofs 

2.1 Dilute Wigner ensemble of random matrices 

Let us consider a family of independent Bernoulli random variables 'Dn,p = 
{dn,p{i,j) ■ 1 < * < i < JT-} with the law 



that is independent of the family of independent random variables An- We 
assume that An and I?„ are defined on the same probability space {fl, F, P) and 
we denote by E{.} the mathematical expectation with respect to P. 
We assume that the random variables a^- satisfy conditions 



where 5ij is the Kronecker symbol. In what follows, we require the existence of 
several more absolute moments of a{i,j) that we denote by 




1 with probability p/n 

with probability 1 — p/n, Q < p <n 



Efly = 0, Ea?- 



(2.1) 



^Xr^ sup E{|a(i, j)!''} 



l<2<j<n 



where the upper bound for r is to be specified. 
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We define the dilute Wigner ensemble as the family of real symmetric n x n 
random matrices Hn,p of the form 

Hn,p{i,j) ^ ^a{i,j)dn,p{i,j), I < i < j < n (2.2) 
VP 

and consider the resolvent 

G„,p(z) = (F„,p-z)-\ Imz^O. (2.3) 

The normalized trace of the resolvent gn,p{z) = Tr Gn,p{z) represents the 
Stieltjes transform of the normalized eigenvalue counting function <t{X; i7„ p) 
(1.2) 

5„,p(z) = - TrG„.p(z) ^ [ ^"^^'^--pI Imz^O. 
n J A — z 

We study asymptotic behavior of gn,p(z) in the limit n,p ^ oo, for z G A^, 

A„ = {zeC: |Imz|>2u + l}. (2.4) 

Our first statement generalizes the result about the semicircle law in dilute 
Wigner ensemble of random matrices obtained under more restrictive conditions 

Theorem 2.1 If the family of random variables An (2-1) is such that ^2+ p < 00 
with p > 0, then gn,p{z) determined by (2.2) and (2.3) converges in probability : 

P- lim 5„ p(z) = w(z), z e (2.5) 

n.p — ^00 

where the function w{z) verifies equation 

wiz) ^ — -— , /rnz^O; (2.6) 

— z — v'^w{z) 

w{z) uniquely determines the semicircle distribution (1.3) being its Stieltjes 
transform and therefore (2.5) implies the weak convergence in probability 

tT(A;iJ„,p) ^ CTsc(A), n,p^oo. 

The proof of Theorem 2.1 is based on the following two asymptotic relations : 
lim E{5„,p(z)} = w(z), z e A„ (2.7) 

n,p — *oo 

and 

Var{(7„.p(z)} = 0(1), z € A^,, as 00. (2-8) 

Indeed, convergence (2.5) can be deduced from (2.7) and (2.8) with the help of 
the standard arguments (see for example [l] or |16)). 

The further improvement of (2.8) is related with the asymptotic properties 
of the resolvent covariance function 

C„,p(zi,Z2) = E{5 (Z2)} - E{5„^p(zi)}E{5„,p(z2)}. 

Let us formulate corresponding statement. 
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Theorem 2.2 Let An be such that, in addition to (2.1), the following properties 
are verified : 

E{a(i,i)3} = E{a(i,jf } = 0, ¥.{a{i,jf}^Vi{l + 5,,)\ (2.9) 
and < oo. Then in the limit — > oo such that 

p^O{n'^), 3/5<a<l, (2.10) 

equality 



C„,p(zi, Z2) = ^5(zi, Z2) + f — - ^ ) T(zi, Z2) + o(n-2) (2.11) 



holds for all zi G Ay with S and T given by the formulas 



1 f Wi- W2 

(1 — v'^w'l){l — v'^w\) \ Zi — Z2 



^(^- = .,2„.,n_..2„.^ ( ) > (2.12) 
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^3^3 

T(zi,Z2) = -- 5-^77; ' (2-13) 

where wi — w{zi) and W2 = w{z2) are the solutions of (2.6). 

Let us discuss results of Theorem 2.2. If one considers the particular case of 
(2.10) when p = n, then (2.12) turns into equality 

C„,p(zi, Z2) = 4 {v^S{zi,Z2) + KiT{z^, Z2)) + o{n-^), (2.14) 

where = V4 — iv^ is the fourth cumulant of the random variable . Relation 
(2.14) coincides with that derived in [TS] for the resolvent covariance function of 
the Wigner ensemble that one gets from Hn,p (2.2) when taking p = n. Therefore 
Theorem 2.2 generalizes the results of |13j . 

It was shown in [TSj that in the local scaling Hmit in the bulk of the spectrum 

zi = A+ + iO, Z2 = A- -iO and Ae(-2w,2w), (2.15) 
2n 2n 

the leading term of (2.14) converges to the expression (cf. (1.7)) 

-| {v^S{zi,Z2) + KiT(z^,Z2)) ^ (2.16) 

where the term with does not contribute. 

In Section 6 we show that the leading term of C„,p(Zi,Z2) (2-11) exhibits 
the same asymptotic behavior in the local scaling hmit (2.15) as the leading 
term of the resolvent covariance function of the Wigner ensemble (2.16). This 
shows that the dilute random matrices considered in the limit of the moderate 
dilution (2.10) belong to the universality class of Wigner (non-diluted) random 
matrices. The lower band 3/5 in (2.10) is due to the technical restrictions related 
with the cumulant expansions we use. We discuss this question in more details 
at the end of the paper. Pushing forward this order of the one could decrease 
the value of the exponent a (2.10). But this demands under more computations 
than of the present paper. 
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2.2 Cumulant expansions and resolvent identities 



We prove Theorem 2.1 and Theorem 2.2 by using the method proposed 
in papers flB] and [16] and further developed in [U |2l [T^ . The basic tools of 
this method are given by the resolvent identities combined with the cumulant 
expansions technique. In present section we present these technical tools and 
explain the scheme of the proofs of Theorems 2.1 and 2.2. 



2.2.1 The cumulant expansions formula 

Let us consider a family {Xt : i = 1, . . . , m} of independent real random 
variables defined on the same probability space such that E{|Xt|'+^} < oo 
for some q G N and t = l,...,m. Then for any complex-valued function 
F{ui, . . . , Urn) of the class Coo(K™) and for all j, one has 



E{X,F(Xi,...,X„)} = ^^E 



[dXtY 



+ egiXt), (2.17) 



where Kr — Cumr{Xt) is the r-th cumulant of and the remainder SqiXt) 
can be estimated by inequality 

f)q+i pirn 

\e,{X,)\ < C, sup I \ > m\X,r'}, (2.18) 

where Cq is a constant. Relations (2.17) and (2.18) can be proved by multiple 
using of the Taylor's formula [T] or by using the characteristic functions method 
(see for example, |13j). 



(X ) 

1) The cumulants Ki- can be expressed in terms of the moments oi X — t : 



Remarks. 
1) The cu 
in particular, 

k[^''> ^ fli, K^^'^ = fl2 - flj, where Mr = E(X[) (2.19) 

Regarding the right-hand side of (2.17) with q — 1, we see that the remainder 
ei{Xt) is given by the following relation 



1 f ,a2^(iJ°M , , 

+ 2^f^^^|- (2-20) 

2) In present paper we are mostly related with the case when E{Xt) — 
E(Xf) = E(Xf ) ^ 0. Then k[^'^ = Kf'^ = K^^'^ = and 

x(^')=/i2, iff *^ = /24 - iff*^ =/i6-15/i4/l2 + 30/li (2.21) 

In this case, the remainders e-q{Xt) of (2.17) considered with g = 1, q = 3 and 
(7 = 5 are as follows : 
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4 

3! 



K. 



{Xt) 
2 

3! 



dXf 



(2.23) 



and 



5! 



--EiX, 



dXf 



dXf 



where for any given i' = 0, . . . ,3, x^.^' is a real random variable that depends on 
Xt and such that < l-'^^fl- In what follows, we omit the superscripts (Xt) 

in the cumulants and use the following denotation 



(2.24) 



dX^ 



dX[, 



2.2.2 Resolvent identities 

Given two n x n matrices A and A such that A~^ and A~^ exist, we have 

A-^ = A-^ - A-^{A- A)A-^ (2.25) 
In the particular case, relation (2.25) leads to the resolvent identity 

{h - ziy'^ = (h - zl^ ^-(h-zl^ ^ (h-h}j{h- zl)~^ (2.26) 

is valid. Regarding (2.26) with h = and denoting G = (h — zl) ^, we get 
equality 



G{i,j) = ^6ij -^J2^{i,s)h{s,j), ^ = -z-\ 



(2.27) 



where h{i,j), i,j = l,...,n are the entries of the matrix h, G{i,j) are the 
entries of the resolvent G and S denotes the Kronecker symbol. 
Using (2.26) we derive for G = {h — zl)~^ , |Im^;| ^ equality 



dG{s,t) _ 1 
dh{j, k) 1 + 6jk 



[G{s,j)G{k,t) + G{s,k)G{j,t)] 



(2.28) 



We will also need two more formulas based on (2.28) ; these are expressions 
for d'^G{i,j)/dh{j,i)'^ and d^G{i,j)/d^h{j,i). We present them later. 
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2.2.3 The scheme of the proof of the semicircle law 

Let us explain the main idea of the proof of Theorem 2.1 that follows the 
lines of the paper p^J. Here we consider the case when ^3 = supj ^ Ejoy ^ < 00. 
Using (2.23) with h — Hn,p and denoting ^ = —z~^, we can write that 



(2.29) 



To compute E{G„^p(i, j)i?,i^p(j, i)}, we use relations (2.17) and (2.28) and 
obtain the following expressions (to simplify formulas we omit here and everyw- 
here below the subscripts n,p when no confusion can arise) ; 

• if j < z 



E{G(*,j)i?(j,*)} -if2(j- *)Em^ 



2 

-E{G(z,jf + G(z,z)G(j-,j)} + e[P 



(2.30) 



with 

.(1) _ 



X2(j,i)E<^iJ(j,i) 



g^G(^j) 



(1)' 



1 



(2.31) 



where we used the denotations of the end of subsection 2.2.1. 
In (2.30), we have used (2.28) in the form 



E 



Also we have taken into account that 

If z < j, then using equality H{j, i) — H{i, j), we get 

dGii,j)\ , (2) 



E{G(z,j)iI(z,j)} = ^2(*,j)E 



dH{t,j) 



V 

n 



E{G{t,jf + G(i,i)G{j,3)} + ef, 
where e\^' is given by (2.31) with Dji replaced by Dij. 



• If j — i, then 



E{G(z, ^)Hi^, ^)} = K, (H(z, z)) E | | + 4f 



2v 

n 



—E{G{^,^n + e^\ 



where e-j is given by (2.31) with Dji replaced by 



(2.32) 



(2.33) 
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Substituting (2.30), (2.32) and (2.33) into (2.29), we obtain equality 

,,2 



E{5} = ^ + ^ E E{G(^, jf + G{i, i)G{j, j)} + e 



(2.34) 



with 



n. ^-^ 



i=l 



J<1 



1<J 



It is not hard to see that the terms ej^ , i = 1,2,3, z G (2-4) are bounded by 
the same variable 2/Z3(|Im zl'^ny^)^^. We present the detailed computations in 
Section 3. Then we can rewrite (2.34) in the form 

E{5n,p(^)} = e + ^v''Y.{gl^^{z)} + ^„,p(z), 

where V'n.p vanishes as n, p — > oo for all z ^ Ky. 

Assuming that the average E{(g(„_p(z)) } factorizes ( see (2.8)), we obtain 
equality 

E{.gn,p(2:)} = C + Ci'^E{.g„,p(z)}2 ^ 

where lim„^p^oo 4'n,p{z) = 0. Then one can conclude that E{(7„^p(z)}) converges 
to the solution of the following equation (cf. (2.6) and (2.7)) : 



w{z) = £,+ ^v'^w{zY . 



2.2.4 The leading terms of resolvent covariance 

In this subsection we present the scheme of the computation of the leading 

terms of C„^p(zi, Z2) (2.11). Let us denote gi = gn_p{z[). I = 1,2. Given a random 
variable, we consider its centered counterpart, /*' = / — E/. Using identity 

E{fg°} = E{fg}, (2.35) 

we rewrite C12 = Cn,p{zi, Z2) as 

n 

C12 = E{5?52} = - E E{5?G2(i, i)}. 

Applying the resolvent identity (2.27) to G2(i,i) = Gn,p{i, j] Z2), we obtain 
equality 

C12 = - ^ E E{5?G2(i,i)if(i,i)}. (2.36) 
To compute E{55'G2(i, j)-ff (j, i)}, we use again (2.17) and get relation 
E{g,G2{^,J)HiJ,^)} = K^E ^ ^^^..^ | + ^ ^ | ^^^..^3 | + 

where ii'r is the r-th cumulant of H{j,i) and Tij vanishes. Using twice (2.28), 
we conclude that 
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dH{j,i) dH{j,z) +^^^''^^n£^ dH{j,z) 

= -g°^[G2{i,jf + G2{i,i)G2{j,i)] - ^Glii,j)G2ii,j). (2.37) 



Then we get equality 



C^2=^2v'E{g',9l} + ^ j2 E{5?G2(i,jr} + ^ ^ E{Gl{i,j)G2{i,j)} 



where Z2), z; G A can be shown to vanish in the Umit n,p — > 00. 

Regarding the right-hand side of (2.38), we apply to the third term the 
resolvent identity (2.25) 

Q Q ^ ^1 ~ G2 

Zl - Z2 

as follows ; 



TrGlG2 = TrG.^^ = ^ - ^^^^^1^. (2.39) 

Zl -Z2 Zl- Z2 {Zl - Z2Y 



Using identity (2.35) that gives equality 

^{glgl) = 2E{ff0g2}E{g2} + ngliglf} (2-40) 
and taking into account (2.39), we rewrite (2.38) in the form 



C12 = 26t;'E{ff2}Ci2 



26^' 



1 E{TrGf } Ejffi} - E{g2} 
n Zl - Z2 {zi - Z2Y 



^4E l^^^-^l Z2). 



(2.41) 



Elementary transformations of the terms in brackets based on convergence (2.7) 
and equality (2.6) allows one to recognize the terms S(zi,Z2) and T(zi,Z2) of 
(2.12) and (2.13) that arise from the second and the third terms of the the 
right-hand side of (2.41). 

3 Proof of Theorem 2.1 

Let us introduce a family of independent real random variables An,'p = 
{«p(*)i) '■ 1 < * < j < defined by 

ap(*,j) = a(*,j)I{K,,,)|<^} = | Q if \a{z,j)\>^, < p < n. 
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where {a{i,j)} verify conditions of Theorem 2.1. We define a real symmetric 
n X n random matrix Hn,p by equality : 

Hn,p{i,j) = -^ap{i,j)dn,p{i,j), 1 < i < j < n 

and consider the resolvent Gn,p{z) = {Hn^p — z)^^ , Imz 7^ 0. 

Regarding gn,p{z) = Ti Gn,p{z), we prove in subsection 3.1 that 

lim Fi{gn,p{z)} = w{z) z e Ay (3.1) 

n,p— ^00 

provided that the variance of gn,p vanishes ; 

Var{,g„^p(z)} = o(l), z e A^, as n,p -^00, z £ Ay. (3.2) 

We prove (3.2) in subsection 3.2. 

At the end of this section we show that 

lim ElgnJz) - gnJz)\ = 0, a.sn,p^oo, z G Ay, (3.3) 

where gn,p{z) is determined by (2.3). Then relations (2.7) and (2.8) follow from 
(3.1) and (3.2) and Theorem 2.1 is proved. 

The proofs of relations (3.1) and (3.2) represent the main subject of this 
section. Let us start to perform this program. The last general remark is that in 
what follows, we will use many times the following two elementary inequalities 

\G{i,j)\ < \\G\\ < (3.4) 

and 

f2\G{i,j)f = \\Ge,\f<-^, i = l,...,n (3.5) 

that hold for the resolvent of any real symmetric matrix. Here and below we 
consider ||e||| = |e(j)P and denote by ||G|| = supngn^^j^ l|G'e||2 the corres- 
ponding operator norm. 



3.1 Main relation for E{^„p(z)} 

Regarding (2.27) with h = Hn,p, we can write that 

E{5n,p(^)} =^--Yl E{G„,p(i,i)^„,p(i,i)}. (3.6) 

To compute 'Ei{Gn,p{i; 3)Hn,p{j., ?')}, wc use formula (2.17) with q = 1 and equa- 
lity (2.28) (everywhere below, we omit the subscripts n,p when no confusion can 
arise). Then we get relation 

E{G(i, j)iJ(j,i)} = k^{j,i)^{G{i,3)} + i^2(j,i)E [DjiG{i,3)} + (3.7) 
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with 



(1) 



DlG{i,j) 



(0) 



where we have denoted Z3Jj = d^/dH{j,iY and Kr{j,i) — Currir (^Hn,p{j,i)j . 

Substituting (3.7) into (3.6) and taking into account formula (2.28), we ob- 
tain equaUty 



z — / 



(3.8) 



where 



R = --Y. ^i(j,»)E{G(z,j)} 



^ [f2(j,») - t;2]E{G(z,jf + G(z,z)G(j,j)} - ^ ^ e,. 



with T/2(j,i) = nK2ij,i)/{l + S,,). 

Now we can rewrite (13.81) in the form 



where 



(3.9) 



E{5} = e + et-'Elsl^ + i? + + 02, 



' -1 



and 



42 = ^v' (E{g2} „ E{g}2 



(3.10) 
(3.11) 
(3.12) 



Let us show that the terms R and 0/, ^ = 1,2, vanish in the limit n,p ^ oo. 

We start with R. Regarding the first term of the right-hand side of p.9p and 
using l|3.4p . one obtains 



J2 Mj,^mG{i,j)}\ 



To estimate the second term of the right-hand side of (|3.9p , we use p.4p , and 
inequality 

\v2{J,^)~v'\<\nai3,^n-nau,^f}\ 



fJ-2+p 



(3.13) 



< E 



|a(j>«)l 



2+p 



P' 



~2 ■^{|a(i.»)l>v^} 



< 



np/2 ■ 



(3.14) 
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Regarding the third term of the right-hand side of (|3.9p and taking into account 
equality 

D%Gii,j) = ^^-i^[2G(z,jf + 6G(z,z)G(j,j)G(*,j)], (3.15) 
and estimate (|3.4p . we conclude that |D^jG(«,j)| < Sjlmzj"^ and that 



,^ E e,.| < 4M ^ fsup|i?|,G(^,J)|) ^imm < (3 



.16) 



Now gathering relations given by (|3.13p . I|3.14p and p.l6p . we get the follo- 
wing bound for R : 

\^\ = o(^) . «s n,p^oo. (3.17) 



Inequality l|3.5p implies that 



|</'i|<^, rj^2v + l. (3.18) 

To estimate 02 (|3.12p . we use the elementary inequality |E{5^} — £{5}^! < 
Var(g„^p(z)) and relation (|3.2p that we prove in the next subsection. 
Relations ([321), (|3T7P and fSTTSl) show that 

\R + (pi + (f)2\ — o{l) , as n,p— >oo. (3.19) 

Then equality (|3.10p and estimate l|3.19p imply that E{5„.p(z)} — > w{z), 
z € Ay, where w{z) is the solution of equation 

w{z) = ^ + ^v'^w{zf , 

such that Imw(z)-Im2 > 0, Imz 7^ 0. This proves convergence in average l|3.6p . 

3.2 Estimate of Var{^„_p(2;)} 

Let us denote gi — n^^ TiGn p{zi), 1 — 1,2. Then we can write relations 
(c.f. (2.36)) 



E{gOgO} = ^dih} = E E{5?G2(z, j)^(j- z)}. 



For each pair (i, j), g\G2(i,j) is a smooth function of i). Its derivatives 
are bounded because of equation (2.28) and l|3.4p . In particular, 



\DlMG2{t,m\ < C (llmzil"! + |Imz2|-^)', 

where G is an absolute constant. 

According to the definition of H and the condition ji^^P < 00 of theorem, 
the third absolute moment of H{j, i) is of order l/{pP/'^n). Then we can apply 
(2.17) with 9 = 1 to ^{glG2{i,j)H{i,i)} and get relation 
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n ^ — ' 

I E i^2(j,.)E{i)), (5^G2(^,J))} - f E 4' (3-20) 



where Kr is the r-th cumulant o{H{i,j) and 



141 < 4G (|Imz,|- + |Imz2r^)^E|i70-,.)P < (3.21) 

Using expression (2.37) and identity (2.40), we rewrite (|3.20p in the form 
E{ff?52} = 26«'E{ff0g2}E{52} +6^^'E{gO(ff°)2} 

. n 

+ ^ E ^2(j,*)E{5?G2(*,jf } + fii2, 

where 

= E ^i(j,0E{g?G2(z,j)}+^ E [1^2(j,*)-«']E{.g?G2(z,z)G2(j,j)} 
Of " t " 



E F2(j,z)E{G?(z,j)G2(*,j)}-f E 4 • (3-22) 
Introducing the auxiUary variable 



n 



we can write the following relation 
E{.g?g2} = g2t''E{.9?(g2")2} + ^ E ^2(j- *)E{g?G2(*, j)'} + 7^^i2. (3.24) 

^ -I S2 

It is easy to see that 

|^2|<^ Z2eA„. (3.25) 



Then 

,.2^7lr;;0/-o^2^ , 92 
n 



|g2«'E{5?(g°)n + ^ E ^2(j, *)E{g?G2(z, j)'}! 

< ^ (Var{gi}Var{.92})'/' + (Var{gi})i/^ (3.26) 
To estimate R12 p.22p . we use inequality 

Y,nG\{^,J)G2{^,])\<[Y.\G\{^,J)\A El^2(^j)n <-5TT 

i=l \i=l / \i=l / ' 

(3.27) 
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with 6 = 2. Then computations similar to those of subsection 3.1 imply that 

1 1 " 



\fRi2\<C 

?2 



pp/2 



(3.28) 



where C is a constant. 

Considering p.24p with z = zi = Z2 and using l|3.26p and (|3.28p . we get 
inequality 



2i]^ 4-1]^ , 

Var{g} < — Var{g} + /Vai^ + C 

77^ nrj-^ 



1 1 



pp/2 j^2 



Then follows. 

Let us prove p.3p . Using the resolvent identity (2.26), we can write that 



,(z)-g„,p(z) = - V G{i,s)[H-H]{s,t)G{t,i)^-y^{GG)(s,t)[H-ms,t) 

11 ^ — ^ n ^ — ^ 



2,s,i — 1 



This relation together with l|3.4p imply that 



1 " 

ngnA^) - 9n.p{z)\ < Y] E|a(s,i)-a(s,t)|E|d(s,t)| 



1 " ]5 f 

E^-E|I«(^'*)|I{KM)I>VP} 



< 



Then relation ll331) follows 



|a(s,i)r+'' I < M2+P 



l+p I — jj2pp/2 ■ 



4 Correlation Function of the Resolvent 

In this section we give the computations that represent the principal part of 
the proof of Theorem 2.2. The auxiliary technical results will be proved in the 
next section. 

4.1 The scheme of the proof of Theorem 2.2 

Let us consider (2.36) and apply (2.17) to E{g^G2{i,j)H{j,i)} with q = 5. 
Taking into account (2.9), we get relation 

^ n ^ n 



120n 



J2 KeI^{D%{g",G2ii,3))}+r, 



(4.1) 



where 
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- I E ^E{i/0-,z)^I?^6j5?G2(*,j)]^"^} . (4.2) 

Let us note that 

,,2 /T/ '^ii'iX A 

if2 = _(l + 5^.,), if4=^_±^ (1+^ = )2 (4.3) 

with 

A = ^4 - 3^4^ 

and 

^6 = 5 — + 1 + ^1'^ = — 1 + -^jO 4.4 

with a = V(i — IdV^v^pn^^ + BOw^p^n^^. In l|4.2p . we have denoted for each pair 



n 



and G'^^^ {zi) = {Hj^ ^ ^ ^i) ^5 ^ = 1, 2 with real symmetric 

^'^'■^^ I if (r-,s) = (j,*) 

where jiJ^'^Hj,*)! < = 0, . . . , 3. 

Regarding the first term of the right-hand side of l|4.ip . we can use (2.37). 
Taking into account (|4.3p . we write that 



Y: E{G?(^,J)G2(^,J)}• 



2J = 1 



Using this equahty and relations (2.39), (2.40) and (4.4), and computing the 
partial derivatives with the help of (2.28), we get the following relation 
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n"' [ (zi - 22) 'T- ^1 - ^2 



26A 



^{Gjii,i)Giij,j)G2{i,i)G2ij,j)} + J2Yr + T + T (4.5) 



r=l 



with 



r2 = ^E(,?iTrG^j, 



= ^ E (E{.,;G2(i.j)*} + liE{9!G2(i,j)'G2(i,i)G2(i,3))), 



^5 = ^ E E{G?(^,J)G2(^,J)' + 3G?(^,J)G2(^,J)G2(^,^)G20•,J)}, 



^6 = ^ E E{G?(*,j)Gi(z,j)G2(*,jf + G?(z,j)Gi(*,j)G2(*,*)G2(j,j)} 



+ ^ E E{G?(z,z)Gi(j,j)G2(*,jr}, 



F7 = 



26A 



E E{G?(*, j)Gi(*, j)'G2 j) + G?(z, z)Gi (j, j)Gi (z, j)G2(*, j)} 



26A 

v?p 



E E{G2(i , j)Gi (j, j)Gi z)G2(i , j)}, 



1^8 = E ( E{5?G2(*, + -E{G?(^, ^)G2(^, i)^ 

^ E f-E{G?(«, *)Gi(z, z)2G2(z, z)} + -E{G?(z, z)Gi(*, z)G2(*, *)^} 



n^p 



and 



T 



6 f:^ii±M!E{i^^^a,0G2(^,,))}, 



120n — ' Tin 



where r is given by l|4.2p . 

Let us discuss the structure of relation l|4.5p . We see that the first term of 
the right-hand side of (|4.5p is expressed in terms of G12. This will finally give 
a closed relation for G12 . The second and third terms of the right-hand side of 
(|4.5p give a non-zero contribution to G12 that provides the expressions of the 
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leading terms <S'(zi, Z2) (2.12) and T(zi, Z2) (2.13). We compute this contribution 
in subsection 4.3. 

The three last terms of (|4.5p contribute as the terms of the order 0(71^^) 
in the Hmit (2.10). The following two statements give the detailed account on 
these vanishing terms. 

Lemma 4.1 Under conditions of Theorem 2.2, the estimates 

\Y,\=0 ({Var(5i)}i/2[p-2 + Var(52)]) , (4.6) 

1^2! - o +p-2{Var(c;i)}i/2 +p-i{Var(5i)}i/2{Var(g2)}i/2) (4.7) 

and 

\Y3\ = O (p-2{Var(gi)}i/2 +p-i{Var(gi)}i/2{Var(g2)}i/2) (4.8) 

are true in the limit n,p 00. 

We postpone the proof of Lemma 4.1 to the next section. 
Lemma 4.2 Under conditions of Theorem 2.2, the estimates 

niax jy^l = O (n- VM""^ + {Var(.gi)}i/2]^ (4.9) 

|T| - O {p-^[n-^ + {Var(5i)}i/2]) (4.10) 
|r| = o (n- V'K' + {Var(.gi)}i/2]) (4.11) 
are true in the limit n,p ^ co. 

Proof of Lemma 4-2. We start with l|4.9p . Inequality l|3.4p and (|3.5p imply 
that if zi e Ay, then 

7A " /I \ 

< E E|5?G2(z, j)'l = O -{Var(5i)}i/2 . 

Inequality ((3^ and fOTP with 6 = 2 imply inequality 1^5 1 < ^cl\l {ifrp-p). 
Using ijOl) . (|3?5l) and l|3.27p with 6 = 1, we obtain that the terms Fe, and 
Is are all of the order indicated in (|4.9p . 

Regarding (|4.10p . it is not hard to see that this result follows from the 
E|i?fa5^G2(^, j)}| = O + {Var(gi)}V2) (4.12) 

in the limit n,p^ 00 (2.10). Let us prove (|4.12p . Using (2.28) and l|3.4p . we get 
relation 

Dji{9°,} = -j2DAGi{t,t)} = --G?(z,j) = O (-) 
n n \n/ 
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for all zi g At,. It is easy to show that 

^;a5?} = O (^i) , r = l,2,..., ze A„. (4.13) 

Then l(4TT2l) follows from l(4TT3l) and (|X4l) . Estimate (|4?T0|) is proved. 

To proceed with estimates of r (|4.2p , then we use the following simple sta- 
tement. 

Lemma 4.3 IJlj If zi e A^, I — 1,2, under condition of Theorem 2.2, the 
estimates 

Var([5„^p(z0]^"^) = O(Var(g„,p(z0)+p-'n-2), = 0, . . . , 3 (4.14) 

and 

l6 



([g°G2(*, j)]('')) = O [n-'\Gi''\t,j)\ + \[9T^Gi''\z, , = 0, . . . , 3 

(4.15) 

are true in the limit n, b — > oo. 

Lemma 4.3 is proved in [I]. We do not present the details here. 

Regarding the first term of the right-hand side of l|4.2p and using l|4.14p and 
(|4.15p . we obtain inequality 



Y^K,miJ,^)DU9'^G,i^,Jm < E S ( ^ + 07^ (Var([5i](^)) 



= o(^ + 4(Var(.gi))'/'), (4.16) 

where c is a constant. Repeating previous computations of l|4.16p . we obtain 
that 

n n 

E^4E|//(j\*)37^,«,[ffOG2(*,jO]^'^l+E^2E|ff(j,z)^I),sj5?G2(z,j)](i)| 

^ + l(Var(.gi))'/'V (4.17) 
Now, regarding the last term of l|4.2p and using (4.15), we obtain inequality 

n 

-EE|H(j-z)'i?fj5?G2(z,j)](°)| 

*,i=i 

< £1 ^ j^E&^ + E|i/(j,z)^[g?](«)||G(")(z,j)|), (4.18) 
i,j=i ^ ' 

where c\ is a constant. Regarding the last term of(4.18) and using l|3.5p and 
(|4.14p . we get 



19 



1 " 1 " I 



4=1 J = l 



1/2 n / n 



1/2 



4=1 \j=l 



.1/2 
Ml4 



<C2^(Var([.9i]("))) ' < ^[(Var(.gi)) 



a/2 



where C2 is a constant. Using this estimate and we rewrite (|4.18p in the form 



-. n . .1/2 / 



a/2 



° — + — (Var(.9i)) 



1/2 



pl/2j^ 



(4.19) 



where c is a constant. Then l|4.1ip follows from the estimates given by relations 
(|4.19p . I|4.16p and (|4.17p . Lemma 4.2 is proved. 

4.2 Estimate of the variance 

Using the definition of q2 p.23p , we rewrite l|4.5p in the form 



Ci 



2^2 



where 



and 



2 ^n,p 



np n 



92 



2 



J2Yr + r + r, 



E{gi}-E{g2} , 1 EjTrG?} 



n zi~ Z2 



= S E E{G?(*,z)Giaj)G2(*,i)G'2(j,j)}. 
Using inequality p.4p and (|3.25p . we obtain that 



414 12^4 

< -g H -g-^, zi,z2 e A^. 



np n 

Lemma 4.1 and Lemma 4.2 together with (|3.25p imply that 



(4.20) 

(4.21) 
(4.22) 

(4.23) 



92 



E^r + T + r 



r=l 



< C 



1 1 

np"^ ri^p 



1 1 

p^ np 



{Var(gi)}V2 



{Var(gi)}i/2Var(g2) + -{Var(gi)}i/2{Var(g2)}^/'') , (4.24) 
P J 
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where c is a constant. Using this inequality and relation (2.10), l|3.27p and i|4.23p . 
we derive form l|4.20p the following estimate 



A 



B 



Var(g„,p(z)) < — W Var(5„,p(z)) + — 

with A and B that depend on z only. Since Var((7„ p(z)) is bounded for all 
2; € A^, then we conclude that 



1 

np 



(4.25) 



Substituting this estimate into (|4.24p . we obtain that 



o 



1 



p^^np) 



This fact together with the restriction (2.10) implies that 



1 1 



p^^Jnp 'n? 



and that the estimate 



^2 ^ 



(4.26) 



holds. This proves (2.11). 

4.3 Leading terms of correlation function 

To obtain the explicit expression for the leading term of C„^p(2:i, Z2), it is 
necessary to study in detail the variables Sn,p and T„.p. Let us formulate the 
corresponding statements. 

Lemma 4.4 If zi e A„, I = 1,2, then under conditions of Theorem 2.2, the 
estimates 



1 



ETrGf 



' 1 - v'^wf 



O 



(4.27) 



1 ^ 3 2 / 1 \ 

- E{G?(z,z)Gi(j,j)G2(*,z)G2(j,j)} = -^^^ + - (4.28) 

hold in the limit n,p — > cx) (2.10). 

Proof of Lemma 4-4- We start with l|4.27p and introduce the variable 

n 

M(z) = Af„.p(z) - - V G{i,j)^. 

n ^ — ^ 
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Applying identity (2.27) to G{i,j) and using formula (2.17) with q = 3, we get 
relation 



E{M(^)} = eE{g(z)} + 2^v^nM{z)g{z)} + ^ G\i, i) + 70 

with 

7o=-^E^4E{i?f,(G2(j-.))} 



v It 

^ 5:E{i7(.,,f [i?f,.(G2(j,.))]<"'} 

c 



n 

J2 K,E{His,j) [Dt^{G'{j,s))f^}, (4.29) 



where Kr are the cumulants of i/(s,j) l|4.3p . Using identity (c.f. (2.35)) 

E{/5} = E{/gO} + E{/}E{5}, (4.30) 
we obtain the following relation for E{Af (z)} : 

E{M(2)} = eE{5(z)} + 2ei;2E{M(z)}E{g(z)} + 70 + 71, (4.31) 

where 

71 = 2ei;2E{Af(z)gO(z)} + ^ ^ G3(z, z). 

1=1 

Relations p.4p and (|4.25p imply the estimate 

Regarding relation l|4.29p and using (2.28) and p.4p . one obtains that 
max ( sup \Dl.G^{s,j)\] < C5 



r=3A 

and that 



- + -l^^':^j\''s,G {s,j)m\H{s,j)\ } < + 

(4.33) 

where C5 is a constant. Relations l|4.32p and (|4.33p imply that 

I70 + 7i| = O ( i ) , as n,p^oo. (4.34) 
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Using this estimate and (2.7), we derive from l|4.3ip relation 

E{M{z)} = Cw{z)[l - 2Cv^w{z)]-^ + O 

Then l|4.27p follows from this relation and equality (2.6). 
Now we prove (|4.28p . Let us consider variable 



1 " 

L(2i,Z2) ^^T. nGli^,^)G^iJ,J)G2{^,^)G2iJ,J)} 
Using l|4.30p . we obtain the following relation for L{z) 

L{zi,Z2) 



i j2 E{G? {^, i)G2 ^)}) E (j, j)G2 (j, j)} j 
1 " 

-EE{G?(^,^)G2(^,^)B?2}, 



n 
1=1 



where 



1 " 

Si2 - -Ee{Gi(j,j)G2(j-,j)}. (4.35) 



n 

To proceed with the estimate of L(zi, Z2), we use the following simple statement 
that we prove in the next section. 

Lemma 4.5 If zi G A.,,, 1 — 1,2, under conditions of Theorem 5.2.2, then the 
estimates 

Var{i3i2} = O ([p-' + (Var{.gi})i/2 + (Varte})^/']') , (4.36) 



and 



E{Bi2} = wiW2 + O (4.37) 

- E E{Gli^, ^)G2{^, ^)} = + O - (4.38) 

hold in the limit n,p ^ co (2.10). 

Now l|4.28p follows from Lemma 4.5 and the definition of L{zi,Z2). Lemma 
4.4 is proved. ■ 

Proof of Theorem 5.2.2. Let us complete the proof of Theorem 2.2. It is easy 
to see that if Z2 G A^, then the definition of 52 l|3.23p . the convergence (2.7) and 
equation (2.6) imply that 

W2 

lim qn,p{z2) = 5— 2 , Z2 G A,,. (4.39) 

Finally, using Lemma 4.4 and relations (2.7) and l|4.39p . we derive from (|4.20p 
relations (2.11), (2.12) and (2.13). Theorem 2.2 is proved. ■ 
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5 Proof of Auxiliary Statement 

The main goal of this section is to prove Lemmas 4.1 and 4.5. 

5.1 Proof of Lemma 4.1 
5.1.1 Estimate of Yi (gini) 

Variable Yi = ^v'^E{g^{g^)'^} (|4.5p admits the obvious bound 

|ri|<^/Vai^yE^. (5.1) 

To proceed with (5.1), we prove the following statement. 

Lemma 5.1 If z <E Ky, then under conditions of Theorem 2.2, the estimate 

E|<?°,p(^)r = O {[p-^ + Var{5„,p(z)}]2) (5.2) 
is true in the limit n,p oo (2.10). 

Let us note that the estimate l|4.6p follows from inequality (|5.ip and l|5.2p . 
Proof of Lemma 5.1. Let us consider the average 

1 " 

with T — 5i5253- We apply to Gi{t,t) the resolvent identity (2.26) and obtain 
relation 

n 

W = -£,aY^ ^{T°Gi{t,s)H{s,t)}. 

Applying (2.17) with q = 2> to E{T"G'4(i, s)i?(s, t)} and taking into account 
(2.28), we get relation 

W ^Uv^nT''{gif} + ^,v'^{^ G4(i,s)2} 
+^ E E Gk{y,s)Gk{t,y)G,{t,s)\ 

{i,3,k) \ t,y,s=l ) 

-|EfE{i?3^(T°G4(M))}+f^ (5.3) 

t,s=l 

with 

^ = E E{i?(s,i)'[i5l(T°G4(i,.))](")}+^ ^ if2E{i7(s,i)3[Z?,^(r"G4(t,.s))](i 

' t,s=l t,s=l 

+ ^ E /^4E{i?(s,t)[7?^,(T°G4(t,5))](^)} , (5.4) 

■ t,s=l 
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where Kr are the cumulants of H{s,t) as in l|4.3p . In l|5.3p . we introduce the 
notation 

^ F{i,j, k) = ^^(1, 2, 3) + 3, 2) + F(2, 3, 1). 

Applying to the first term of the RHS of l|5.3p relation (2.39) and using the 
definition of l|3.23p , we obtain that 

^0 " 



where 



t,y,s=l 



^ = --E^E{Z?3^(T°G4(M))} + |if^. 

n ^-^ 6 ?4 



(5.5) 



Regarding Gk{t, ■) and G4(i, •) in the third term of the RHS of (??) as a vectors 
in n-dimensional space, we derive from estimate (|3.5p that 



E Gk{y.s)Gu{t,y)Gi{t,s) 
y,s=i 

<\\Gk\\(J2\Gk{t,y)\'] [Y.\G,{t,s)\' 

\y=l ) \s=\ 



< 



(5.6) 



Now gathering relation given by ijOl) . (|3?5|) . l|3.25p . and 

E|T"(g°)2|<-[E|T|E|g2l+E|r52|] 
imply the following inequality for W : 

47)2 A 2 A 2 1 o 2 

IVFI < — E|r54"| + — E|T|E|5°| + — E|r| + ^E\g^g",\ + (5.7) 

Henceforth, for sake of clarity, we consider G = Gi G3 = G4 = G2, then we 
get T = (5°)2g° and 



Let us assume for the moment that 

1 1 ^Var(.g) Var(.9) ^/Yar{g)VW 



(5.8) 



np 



(5.9) 



Now returning to (|5.7p and gathering estimates given by relations l|3.4p , (|5.8I 
and l|5.9p imply the following estimate 

W <aJ- + v/Var(5)) + — f i + /vMff) 



np \p 
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where Ai, A2 are a constants. Then we obtain l|5.2p . 

To complete the proof of Lemma 5.1, let us prove l|5.9p . To do this, we use 
the following statements. 

Lemma 5.2 If z <E A^, then under conditions of Theorem 2.2 the estimates 

DUT'Git, s)} = O + n-2|50| + n-^lg^ + \T°\) , (5.10) 

[Dt,{T''G{t,s))]^''^ =0[{n-'+n-'\[gT^\+n-^\[g"](^^ 

(5.11) 

and 

E|[gO]Mr = o(p-'-/2n-'- + E|5°r), r-l,...,4 (5.12) 
are true in the limit n,b ^ 00 satisfying (2.10) and for all = 0, 1, 2. 

We prove this Lemma at the end of this subsection. 

Let us return to the proof of l|5.9p . Regarding the first term of the RHS of 
(|5.5p and using the definition of l|4.3p , inequality l|3.25p and estimate l|5.10p , 
one gets with the help of (|5.8p that 

" 9 A 

l|E^Mfi(T»<!(M))}l 



^ 1 ^7Var(i)^Va^^vA^VW 

\ pn^ pn^ np p I 



Now let us estimate O l|5.4p . Regarding the first term of the RHS of (|5.4p and 
using l|5.1ip . we obtain inequality 



1 " 



n 

t,s=l 



< 



n \ n'^ n? n ' 



^^E|i7(s,i)5||[g°](°)nGW(t,.) 



To estimate the last term of this inequality, we use l|3.5p and l|5.12p , and we get 
estimate 

1 " 1 " ^ / 

n — ' n — ' » 

t,s=l t,s=l 
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Using this estimate, relation (|5.12p and the same arguments in the proof of 
estimate of l|4.1ip . one obtains that 



1 " 

- J2 m{s,tnDt,{T''G{s,m 



(0)1 



n 

t.s=l 



-0(^.4^,^.^.^^^Sm. ,5.4) 

Repeating the arguments used to prove l|5.14p . it is easy to show that the term 
1 " 

- J2 i^4E|77(,s,i)[Z?4(T0G(s,t))](i)|+i^2E|i7(s,i)3[i?4^(r°G(s,i))](2)| 
n ^ — ' 

t,s=l 

is of the order indicated in the RHS in (|5.14p and that 



Then the estimate l|5.9p follows from l|5.13p and l|5.15p . Lemma 5.1 is proved. 

Proof of Lemma 5.2. We start with (|5.10p . Remembering that T — [g°]'^g° 
and using (2.28) and (|4.13p . we obtain that 

B\Dl,{T'}\^0{n-'\gr)- 
Dim = 0{n-'\9'\+n-'\gr), 

DUT"}^0{n-' + n-^\g°\+n-'\gr) 
Now it is easy to show that l|5.10p is true. 

Similar computations prove the estimate IjS.lip . 

Finally, let us prove (|5.12p . To simplify computation, we denote [g]^"^ — g. 
Then the resolvent identity (2.26) implies that 

^=-E^(^'^) = -E^(^'^)-- E G{k,r){H-H}{r,^)G{^,k) 

k—1 k—1 fc,r,i— 1 

= g--Tv{GG6H) 
n 



with 

SHir,i)={H-H}{r,t) = 
where < \H{s,j)\ < \H{s,j)\. Then 



if (r,i)^(s,j) 

H{s,])-H{s,j) if (r,i) = (s,j), 



< cE|.g°r + —E\Tr{GG6H) - I^{Tr{GG5HW, r = 1, . . . ,4, (5.16) 

n' 
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where c is a constant. Using l|3.4p . we obtain that 

E\TriGGS„) - E{Tr{GGS„))\ < ^E{\H{s,j)\ + E(|H(s, j)|)} 



(5.17) 

Relation l|5.12p follows from l|5.16p and estimate l|5.17p . Lemma 5.2 is proved. 



5.1.2 Estimate of ^2 (mil) 

Remembering that Y2 = ^2v'^n~^ s ^{diTfGl}, we consider the average 



Y2 = \Y. nglG^ii.sf) 

z,s — 1 

and apply to G2{i,s) the resolvent identity (2.26). Then 

Y2 = ^Cu~^ V E{9lG2{^,s)G2{^,t)Hit,s)}. 

71 ^ — ^ 

i,s,t—l 

Applying (2.17) with g = 1 to 'E{g'(G2{i, 3)G2{i,t)H{t, s)}, we obtain that 



Y2 = 2(2v'E{g2}Y2 + |Ci2 + ^ E{g?G2(*, s)^^^} + ^ e., (5.18) 



i,s — l r— 1 

where 

01 = ^ E E{G?(.,t)G^(.,t)} + ^^$:E{3°G3(.,.)}, 

s,t=l s=l 

©^--|EfE{i^MG^(M))} 

t,s=l 



and 



' t,s=l 

|5:X2E{i7(M)^[^t(3?G^(i,^))](^^} 



+ E ^4E{H(t,.)[i?f,(5?G^(t,.))](2)}, 

■ t,s=l 

where i^r are the cumulants of H{t, s) as in l|4.3p . 

The term 2^2V^'E{g2}Y2 can be put to the left-hand side of (|5.18p . Using 
(|3.4p . I|3.5p and l|3.27p . it is easy to show that the second and the third terms 



28 



of the RHS of l|5.18p and 61 are of the order indicated in the RHS of l(4?7l) . 
Using similar arguments as those of the proof of l|4.6p (see (4.16)-(4.19)) and 
the following estimate (c.f. (|4.12p ) 

DMGl(t,s))^0{n-^ + \g',\), r = 3,4, 

we conclude that the terms 62 and 63 are of the order indicated in the RHS of 
(liJj) . Relation l(4?7l) is proved. ■ 

5.1.3 Estimate of Fa (mH) 

We rewrite I3 in the form Y3 = C2Ap-iE {5?[B22]^}, where (cf. (035])) 

n 

B22 = - VG2(i,z)^ 

n ^ — ' 

1=1 

Let us note that the estimate l|4.8p follows from (2.40), inequality 

|E{5?(i322)'}| < 2E|g°i?°2|E|S22| +E|.g?(i?2°2)'l < -4/^^aiTOv/Var{i?22} 
and estimate (|4.36p . This proves l|4.8p . Lemma 4.1 is proved. 

■ 

5.2 Proof of Lemma 4.5. 

5.2.1 Estimate of Var{Bi2} dUM]) 

Let us consider the average 11 = E{i3]'2G'i(j, j)G2(j, j)} and apply 

to G2{j,j) the resolvent identity (2.26), we obtain that 

n = C2E{B?25i} - ^ E E{i3?2Gi(j,j)G2(j,s)iJ(s,j)}- 

Now applying formulas (2.17) with q — 3 and taking into account relation (2.28), 
we obtain that 

n =C2E{i3025i} + C2^;'nE{g2} + C2V^-E{B°^Buf2} 

7 2" 

+ ^ E E{S?2[Gl(j,j)G2(j,5)2+2Gi(j,j)Gi(j,5)G2(j,5)]} 
972 " 

+ ^ E E{Gl(^,^)G2(^,s)G2(^,J)Gl(J■ j)G2(j,5)} 

+ ^ E E{Gi(*,5)Gi(z,j)G2(*,0Gi(j- j)G2(j,s)} + C/(z) (5.19) 
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with 



j,s=l 



7 

-|j ^4E{ff(.,j)[I?^^.(i3?2Gi(j,j)G20-,5))](')} (5.20) 



2,4 are the cumulant of H{s,j) as in (|4.3p . Let us assume for 
the moment that 



where Kr, r 



U{z) 



o(^ + JVVar{Bi,}). 



(5.21) 



Now returning to (|5.19p and gathering estimates given by relations p.4p , l|3.5p , 
(|5.6p and (|5.2ip imply the following estimate 



Var{Bi2} < Ai 



v/Var{Bi2} 



where Ax and are some constants. Then (|4.36p follows from this inequality. 

Now, let us prove l|5.2ip . To do this we use the following statement. 
Lemma 5.3 If z E A^, then under conditions of Theorem 5.2.2, the estimates 



DlJB'>,Gx{j,j)G2{j,s)}^0i 



Var{[Bi2]('''} = O {p-\i-^ + Var{Bi2}) 
are true in the limit n,p ^ co satisfying (2.10) and for all = 0, 1, 2. 



(5.22) 
(5.23) 



We prove this Lemma at the end of this subsection. 

Let us return to the proof of l|5.2ip . Using the definition l|4.3p and esti- 
mate l|5.22p . it is easy to show that the first term of the RHS of (|5.20p is of the 
order indicated in the RHS of (|5.2ip . 

Regarding the first term of the RHS of (lOO]) and using l(532l) and l(533l) . 
we obtain inequality 



1 " I 

- ^ ^\H[s,jf[D%{Bl^G,{j,j)G2[j,s))] 



(0) 



< 



c 



n 



EE 



n 



Ai5 



' 1/2 

c Mio 



1 



= O 



j-s=l 

1 1 
— + - 



pl/2^ 



v/Var{Bi2} 



VVar{Bi2} , 



(5.24) 
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where c and c are some constants. Repeating the arguments used to prove 
(|5.24p . it is easy to show that the third and the fourth terms of the RHS of 
([5?20l) are of the order indicated in the RHS of (fOTj) . This proves fOTI) . 

Proof of Lemma 5.3. The estimate l(532l) follows from (2.28), ([331), ([33]) and 
(13:271) . 

Let us prove l|5.23p . To simplify computation, we denote [G;]*^"^^ — Gi, I — 
1,2. Then the resolvent identity (2.26) imply that 



Gi{k,k)^Gi{k,k)~ Giik,r)6Hir,i)Giii,k) 

= Gi{k,k)- Gi (fc, s) [H{s, j) - H{s, j)]Gi (j, k) 

for ^ = 1, 2 with 

'"^'^'^-^''-''^^'^''^-[His,j)-His,j) if (r,.) = (.,,), 
where < \H{s,j)\ < \H{s,j)\. Then 
1 " 

Bi2 = - y2Giik,k)G2ik,k) 

k=l 

1 " 

= B12 - - V Gi(fc, s)G2(fc, s)Gi{j, fc)G2(j, fc)<5//(s, 
n -"^ — ' 

fe=i 

1 " 

- - V Gi(A;, fc)G2(A;, s)G2(j, fc)5H(s, j) 
1 " 

- V G2{k, k)Gi{k, s)Gi(j, fc)Jff (s, j). 

rj ^ — ^ 
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fc=l 



It is not hard to see that the last equality together with relations (|3.4p and 
((3:271) implies that 

Var{Bi2} < 4Var{Bi2} + ^ (E|i/(s, j)!^ + [E|i/(s, j)p]% 2E|i/(s, j)p) 

<4Var{i?,2} + 4f^ + 4 + — 

where where ci and c are constants. This proves (|5.23p . Lemma 5.3 is proved 
and this proves of the estimate (|4.36p . ■ 

5.2.2 Proof of relation (gUll) 

Remembering that B12 = rt^^ Gi (i, i)G2(i, ?). Applying relation (2.26) 
to one of G2(j, i) and using formula (2.17) with q = 3, we get relation 

E{Si2} = C2E{gi} + C2w'E{Bi252} + Ti + r2 
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with 

C2V 



Fi =^ V I^{Gi{i,s)G2ii,s)G2{i,i) + Gi{i,i)G2{i,s)^} 



i,s—l 



C 

^ ^ X4E{i?3^(Gl(^,^)G2(^,s))} 

-i.s— 1 



^ V f- 



i,s — 1 



where 



|f..| <csup|Df,(Gi(*,z)G2(*,.s))|E|iJ(s,z)|5 < "^^^ 



iCr- are the cumulants of H{s, i) as in l|4.3p and c is a constant. 

Using identity l|4.30p . we obtain the following relation for E{Bi2} : 

E{Bi2} = ^ , + ^ , [C2^''E{Bi2.92"} + Ti + r2] (5.25) 

Relations p.4p . I|4.25p and the estimate 



max sup|L»^,(Gi(i,i)G2(i,s))| < ci 
'-=3,4 V s,i J 

imply that 

|C2«'E{Bi2.g^} + Ti + r2 1 = O Q j , (5.26) 

where ci is a constant. 

Now to proceed with the estimate of the first term of the RHS of (|5.25p . we 
use the following simple statement that we prove in the end of this subsection. 

Lemma 5.4 If z £ A^, under conditions of Theorem 2.2, then the estimate 

E{5n,p(2)} = w{z) + O (5.27) 

holds for enough n, p satisfying (2.10). 

Using Lemma 5.4 and relation l|5.26p . we derive from l|5.25p estimate l|4.37p . 

5.2.3 Proof of relation (gllHl) 

We introduce the variable 



n 

Ui2 = -y^Gl{i,i)G2{i,i). 

n. ^ — ^ 



n 

i=l 

Applying relation (2.26) to one of G2(i, i) and using formula (2.17) with g = 3, 
we get relation 

E{i7i2} = C2-E{Tr Gj} + C2V^'E{Ui2g2} + *i + *2 
n 
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with 

,,2 



i.s — 1 
m2 



n ^ — ^ 

where 



Z.S — 1 — 1 

n ^ — ' 

1 



|*,.| < csup \D%XG\(^.^)G2{^, s))|E|77(s, < 



ifr are the cumulants of H{s, i) as in l|4.3p and c is a constant. 

Using identity l|4.30p . we obtain the following relation for E{Ui2} : 

Relations p.4p . I|4.25p and the estimate 

anp\Dr{Gl{i,i)G2{i,s))\) < c. 



max 

r=3,4 



imply that 

\C2V^E{Ui292} + Ti + - O , (5.29) 

where ci is a constant. 

Using relations (|4.27p . (|5.29p and l|5.27p . we derive from (|5.28p the estimate 
(|4.38p . Lemma 4.5 is proved. 

Proof of Lemma 5.4 

Remembering that g — n^^J2iG{hi)- Applying relation (2.26) to one of 
G{i,i) and using formula (2.17) with g = 3, we get relation 



with 



E{5} = C + C«'E{g}2 + $i + $2 (5.30) 



i,s — l 

C 



6„ E^4E{Z)^,(G(z,.))} 

$2 = - - V 

71 ^ ^ 



i,s—l 

n ^ 

i,s—l 



where 



|$..| <csup|I?^^(G(*,s))|E|iJ(s,z)|5 < 
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Kr are the cumulants of H{s,i) as in l|4.3p and c is a constant. 
Relations p.4p , (|4.25p and the estimate 



max s\xp\D''^^{G{i,s))\ < a 
'-=3,4 \ s,t J 

imply that 

|*i + *2|=oQ), (5.31) 

where ci is a constant. Using relation (|5.3ip . we derive from l|5.30p the esti- 
mate (|5.27p . Lemma 5.4 is proved. ■ 



6 Scaling Limit And Universality Conjecture 

The asymptotic expression for C„^p(zi,Z2) obtained in Theorem 2.2 and 
T{zi, z-i) regarded in the limit zi = Ai + iei, Z2 = A2 + i(-2 and | 0, j = 1, 2 
can supply the information about the local properties of eigenvalue distribution. 
We follow the schema proposed in [13] . 

Let us recall the inversion formula of the Stieltjes transform (2.6) of 

the semicircle distribution with the density psc = : 

Psc(A) ~ TT^^ limImw(A + ie) = Ix{w{z)}. (6-1) 

eiO 

Consider the density-density correlation function 

r„,p(Ai, A2) - E{p„.p(Ai)p„,p(A2)} - E{p„,p(Ai)}E{p„,p(A2)}. (6.2) 
It is easy to see that the Stieltjes transform of F„,p (Ai, A2) is 

C„,p(zi, Z2) [ [ ^".p(-^i'-^2) — ^XidX2, Imz^T^O. 
J J (Ai - zi)(A2 - 2:2) 

Applying formally the inversion formula (|6.ip . we obtain the following relation 

r„,p(Ai,A2) = /ai O Ix^{Cn,p{zi,Z2)}. (6.3) 

In Theorem 2.2, we have found explicitly the leading term of Cn,p{zi, 22) in the 
domain | Imzj > 2v. However, since the functions S (2.12) and T (2.13) can be 
continued up to the real axis with respect to the both variables zi and 2:2 , we can 
apply to the leading term of (2.11) the operation Ai 7^ A2 to compute 

formally the "leading" term of the density-density correlation function. This 
means that we perform first the Hmit n,p ^ 00 and then the Hmits ei,e2 i 0. 
This order of limiting transitions is inverse with respect to that prescribed by 
the definition (|6.3p . 

Let us denote wj — w{zj), j ^ 1,2 and write the identity 

Wi — W2 'W1W2 
21—22 1 — V^WiW2 

that is an easy consequence of the equation (2.6). This identity yields relations 
e|w(A + ze)p = Imw(A -I- ie){l - v^\w{X + ie)^) (6.4) 
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and |w(A + ie)p — v for A such that Im w(A + iO) > 0. Combining i|6.4p with 
(1.3), we obtain that 

\2 \2 

v^\Rew{\ + iO)Y = — and v^[\mw{\ + = I - —. (6.5) 

Let us consider the terms of the RHS of (2.11) given by (2.12) and (2.13). Using 
(|6.5p . we obtain that 

. . r^l X, 1 4t;^ - A1A2 

iAi o iA.iD^zi, Z2;| 2^;%2[(Ai - Az)]^ (4z;2 - \ly/^{4v^ - X^y/^ 

and 

T T iT( 1 (2z;^-A^)(2t;^-Ai) 

^ 47r2«8 (4i;2 - A2)i/2(4i;2 _ A2)i/2 ' 

Then for Fn,p(Ai, A2) (|6.2p we get the following formal expression 

1 4w2_AiA2 



Fn,p(Al, A2) 



7r2[n(Ai - A2)]2 (4^2 - \lY/^{Av^ - A2)i/2 
1 Vi (2w2-A?)(2w2- A2) 1 3^4 (2i;2- A?)(2w2- A2) 



np 27r2w8 (4u2 - A2)i/2(4„2 _ a2)i/2 n2 27r2w8 (4t>2 - A2)i/2(4„2 _ ^2)1/2 • 

(6.6) 

It is easy to see that in the scaling limit 

Ai,A2->A, n(A2-Ai)^s, (6.7) 

one gets equality 

lim r„,j,(Ai, A2) = 5-^. (6.8) 

We see that the terms of the order 0{l/np) that depend on the value V4 disap- 
pear in the local scale limit (6.7) and the density-density correlation function 
gets the universal form (6.8). This result can be regarded as an evidence of the 
fact that the moderate dilution of Wigner random matrices given by the limit 
(2.10) does not change the universality class of the local eigenvalue statistics. 

Condition 3/5 < a < 1 (2.10) of Theorem 2.2 is related with the technical 
restriction of the approach we use. Indeed, we need this in the bound (4.26) to 
estimate the term T (4.10) that corresponds to the last term of the cumulant 
expansion of the order 5 (see (4.5)). Pushing forward this expansion to the orders 
higher than 5, one could consider lower values for a. However, this requires much 
more cumbersome computations than those of the present paper. 
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